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Abstract 

In this paper we continue the study of the truncated conformal space approach 
to perturbed conformal field theories, this time applied to bulk perturbations 
and focusing on the leading truncation-dependent corrections to the spectrum. 
We find expressions for the leading terms in the ground state energy divergence, 
the coupling constant renormalisation and the energy rescaling. We apply these 
methods to problems treated in two seminal papers and show how these RG 
improvements greatly increase the predictive power of the TCSA approach. One 
important outcome is that the TCSA spectrum of excitations is predicted not to 
converge for perturbations of conformal weight greater than 3/4, but the ratios of 
excitation energies should converge. 
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1 Introduction 



The Truncated Conformal Space approach (TCSA) of Yurov and Zamolodchikov Q has been 
a widely-used method to study the finite-size dependence of perturbed two-dimensional con- 
formal field theories for quite some time ^, ||, |^. It is based on truncating the 
infinite dimensional Hilbert space to a finite-dimensional system on which the Hamiltonian 
is studied numerically-Qlt has been known for a long time that the method has various con- 
vergence problems which can reduce its effectiveness |2|, |3| . The principal problems that have 
been noted before are the divergence of, and differences between ground state contributions 
in different sectors. Once these have been taken into account by considering only differences 
of energy levels in the same sector, the effects of truncation can still be important. There has 
been some interesting work on extrapolation in truncation level 0, |5| , assuming an underly- 
ing but unknown scaling behaviour. Here, we conjecture that the most important corrections 
after the ground state divergence are a renormalisation of the coupling constant and a renor- 
malisation of the energy scale. Our main results are perturbative expressions for the leading 
coupling-constant and energy-scale renormalisations. We apply these to the tri-critical Ising 
model as considered in both and and show that the behaviour of the TCSA results is 
greatly improved. We also find numerical results for the renormalisation and rescaling for the 
tri-critical Ising model and show that the TCSA approach remains very accurate even when 
these two effects become large. 

One important result we find is that while the TCSA estimates of energy gaps converge 
for perturbations with conformal weight less than 3/4, they do not for those with weight 
greater than 3/4. We show this in the case of the minimal model Mg^io perturbed by a field 
of weight 4/5. This failure of convergence is entirely due to the divergence of the energy 
rescaling; once the energy rescaling is taken into account the TCSA estimates converge as 
well as before. In other words, the TCSA estimates of ratios of energy gaps converge even if 
the gaps themselves do not. 

The paper is organised as follows. We first introduce the TCSA approximation and present 
the problems to be addressed in section |2|. We then discuss the perturbative results for the 
coupling-constant and energy renormalisations in sections ^ and ^, and present the results for 
the models considered in |||] and in sections 0| and In section ^ we then discuss the 
problems occurring with /i > 3/4 and in section the complications for perturbed models on 
the strip. 



2 The TCSA approach and its errors 
2.1 The TCSA approach to bulk perturbations 

At present, the TCSA can be applied to bulk perturbations in two different arenas: when the 
model is defined on a circle and when it is defined on a strip. The original paper of Yurov 



and Zamolodchikov treats the cylinder case, the boundary case was initiated in [11]. We 
consider here the case of the cylinder, and defer the strip to section |7[ 

We start with a CFT defined on a cylinder of circumference R which we take to be a strip 
< y < i? of width R in the complex plane with coordinate z = x + iy, and with the two 
edges of this strip identified. For the details of CFT, see [12|. 



^In this paper we study the original form due to Yurov and Zamolodchikov, not the revised version of 
0. 
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The unperturbed Hamiltonian generating translations along the cylinder is 



H= Tr,, ^ . (2.1) 

Jo 

We will map the strip to the complex plane with coordinate w = exp(27rz/i?) in terms of 
which the CFT Hamiltonian is 

H=^[Lo + Lo-^) , (2.2) 

where Lq and Lq are the zero modes of the two copies of the Virasoro algebra present in the 
theory defined on a plane. 

We are interested in perturbations by one or more bulk fields (pi{x,y). We take these to 
be spinless quasi-primary fields of equal left and right conformal dimensions {hi, hi). If the 
coupling to these fields are /Uj then the perturbation is given by an addition to the action 

5S= f fornix) d^^- (2-3) 

i 

When mapped to the upper half plane this gives the perturbation to the Hamiltonian as 

6H = Y,^^^[^) J^_^V^{e'')de, (2.4) 

where w = rexp{i9) so that y = {2-k9/R). 

Since the circle has rotational symmetry, one can restrict attention to the rotationally 
invariant states on which Lq — Lq = 0. On these states, we can perform the ^-integral in 
(2^) so that the perturbed Hamiltonian becomes 

2tt 



^ R 



Lq + Lq-^ + J^(27r)^^i'' a, 99,(1). 



(2.5) 



where = 2 — 2/ij and Aj = ^iR^"- are dimensionless coupling constants. This is the TCSA 
Hamiltonian of Q. 

The TCSA method is to truncate the Hilbert space in some prescribed manner and take 
the Hamiltonian to be ( |2.5| ) on the truncated space. There are at least two ways this has 
been implemented. 



2.1.1 Truncation by level 

The first method is to truncate the space to level n in each representation in the Hilbert 
space, as in [|l|. This means that the maximum value of Lq + Lq in each representation is 
(2n + hi + hi). This can cause problems if the Hilbert space includes some representations 
with large values of /i, as these high level states cause distinct qualitative changes to the 
spectrum. 
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2.1.2 Truncation by total energy 

The second method is to truncate so that Lq + Lq < 2n; this means that the maximum 
value of Lq + Lq in each representation is 2n, and so if a representation has /i > n, it will be 
completed excluded from the TCSA space at level n. This method has the advantage that 
high weight representations do not unduly affect the TCSA space, that the space at level n is 
smaller than in the truncation by level, but has the disadvantage that the truncation affects 
are harder to deal with analytically. 

These two methods will give very similar results when n ^ h for all /i in a particular model. 
It is worth pointing out that choosing a different truncation method can have dramatic effects 
- the "mode truncation" investigated by Toth in has a very different behaviour. We shall 
use the level truncation unless stated otherwise. 

The operator H in ( p.5| ) is dimensionful, and it is far preferable to work with dimensionless 
operators, or dimensionless eigenvalues. In the case of bulk perturbations, there are two 
natural choices. 



2.1.3 Flovi^s ending in massless theories 

If the IR limit of the flow is a massless theory, then it is natural to work with the scaling 
functions defined in terms of the energy eigenvalues £"„ as 

en{^^) = (^] En{K) , (2.6) 



2tt^ 

The scaling functions are expected to flow to the eigenvalues of the operator 

(Lo + Zo-^),^ , (2.7) 
at the IR fixed point with a correspondingly simple spectrum. 



2.1.4 flows ending in massive theories 

If the IR limit of the flow is a massive theory, then the mass of the lightest stable particle, 
m, gives a natural scale and one considers the dimensionless operator 



H _2-K 
m r 



A,: 



(2. 



where r = mR is a dimensionless variable and A,m are a set of numbers. 



2.1.5 Problems vi^ith the TCSA 

The original model studied with the TCSA by Yurov and Zamolodchikov was the Lee- Yang 
model perturbed by its primary fleld v?(i,3) of weight h = —1/5. The results in this model, 
whether as in the original case on the cylinder or on the strip |jll|, have been exceptionally 
accurate, up to 14 digits for some quantities. When the TCSA method was applied more 
generally in and Q, it was clear that this was not always the case. There were several 
problems identified in Q for the perturbation of the tri-critical Ising model by the primary 
field of weight 3/5; we discuss these problems in turn. 
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The first is that the Hilbert space of the model may split up into several sectors and the 
TCSA eigenvalues in these sectors may differ by unphysical amounts. In our prototypical 
examples of perturbations by Lpi^, and ip^i, the perturbation will only couple together repre- 
sentations in the same row or same column of the Kac table. For example, the tri-critical 
Ising model has representations {(r, s); 1 < r < 3, 1 < s < 4, r + s even}. Under the action of 
(^(13) these fall into three sectors which we denote (r, *) for r = 1, 2 and 3. Under the action 
of ¥'(3,1); they split into four sectors, which we denote for r = 1,2,3 and 4. Under the 

massless flow generated by ^(1^3), the (*,r) sector flows into the (l,r) representation of the 
Ising model. Under the massive flow generated by <P{i^3), the three sectors flow into three 
sectors of the massive model, corresponding to the splitting of the three degenerate ground 
states of the massive kink model. Under the irrelevant flow generated by V'(3,i)) the (*,r) 
sector flows into the corresponding (r, *) sector of the minimal model M(5,6). The problem 
is that the TCSA gives slightly different results for the ground state energy is the different 
sectors, so that it may be impossible to get accurate results for energy differences of states in 
different sectors. Sometimes interpolation in truncation level or size of the truncated space 
produces sensible results, as was reported in [18|. The main result we have to report is that 
the energy rescaling formula we find in section ^ is slightly different in different representa- 
tions, and that this effect is most noticeable for the ground states in each sector. This greatly 
improves the difference between the ground states in the different sectors, as can be seen in 
figure |5[ 

Secondly, if the weight of the perturbing field becomes larger than 1/2, the conformal 
perturbation expansion and correspondingly the TCSA eigenvalues, become divergent. This 
has the result that the ground state energy of the TCSA system does not converge with 
increasing truncation level. This was observed in and discussed in [^]. It was realised 
that this divergence can arise purely from the second order contribution for which there is an 
exact expression. Subtracting the divergent part of this expression then gives a revised TCSA 
estimate which will converge, with increasing truncation level, to the perturbed conformal 
field theory result, as observed by Takacs [^l|. We re-derive this leading term as part of 
our treatment of the energy rescaling in section ^. We illustrate the effectiveness of these 
subtractions in the case of the massless perturbation of the tri-critical Ising model considered 
in 0. We show in figure 1(b) that this works well in the case of the tri-critical Ising model 
perturbation considered problematic in - after subtraction of the leading divergence, the 
ground state energy does then converge for this perturbation. 

Finally, it also appears that the "scaling region" is not easily reached, the region where 
the eigenvalues scale with truncation level in the expected manner. Our solution is that 
suggested in a careful consideration of the scaling of the model with system size which 
we show reduces, in the cases considered in [Q], to a renormalisation of the coupling constant 
and a representation-dependent re-scaling of the Hamiltonian. This is almost the same as we 
found in [16 1 . As in |jl^, perturbative expansions for these effects can be found by considering 
the change in the energy eigenvalues with truncation level. We consider first the coupling 
constant renormalisation in section ^, and then the ground state energy and energy rescaling 
in section ^ and apply them to the tri-critical Ising model in section ||. 
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3 Coupling constant renormalisation 

The derivation of the couphng constant renormahsation is a straightforward generalisation of 
the boundary case. We assume that the perturbed correlation functions on the cylinder with 
coordinates (x, y) are given by the insertion of the expression 



(i-oo r-R 
-^f^i / Pnfi{x, y)cyl Pn dyd 



xj , (3.1) 

in the unperturbed expressions, where V denotes path ordering and P„ denotes the projector 
onto states at level n or lower. After mapping to the plane with w = exp{2Trz/R), r = 
exp{2'Kx/R), 6 = 2Try/R, this becomes 



(poo /'2n 
- VAi(27r)-J'' / / Pn 
Jr=o Je=o 



rdrdO 



= / Pn^i{w,w)cyl. Pn -—r + ■ ■ ■ , (3.2) 

Jr=o Je=o 

where yi = 2 — 2/ij and Aj = (27r)~^'Aj is introduced for convenience. We require that the 
perturbed correlation functions be invariant when the truncation level n is changed. The 
simplest way to find the leading order change in the coupling constants is to consider the 
matrix elements of the integrand of r^~^dr in ( |3.2| ) taken at r = 1 and taken in the states 
{ifil . . . |0), that is we consider 

^ ^ fi r2iT drdO 

Zi,„ = -(27r)Ai + ^(2^)A,Afc / / ((^,|^,(1, l)P„(^fc(r, 0)|O) --^ + . . . , (3.3) 
j^f^ Jr=o Je=o 

where we have performed one of the angular integrations. We have also used the first of the 
properties of the primary fields 



{ipi\(fj{w,w)\0) = 6ij\w\ 



-2K 



(^,|V.,(1,1)¥^.(^,U))|0) = (3.4) 



Requiring = we find 

'1 r2-n 



f f drd^ 
A,(n)-A,(n-1) = VAjAfe / / {ipi\ipj{l,l)[Pn - Pn~i]Mr,0)\O)—ZT ■ (3-5) 
^ Jr=o Je=o 

From (U), 

T{hj + hk -hi + n) 







1) [Pn - Pn-l] ^k{r,e) \^)dd = 27rC,jfc 



T{hj + hk- hi)T{n + 1) 



r^". (3.6) 



Substituting (3.6) in (|3.5|), performing the r integral and expanding out to leading order in 
n, we get 

jk 
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As we see, there are corrections to Aj from all pairs of fields <y9j, ipk which couple to ipi, but 
that those for which hi — hj — > 0, i.e. those which appear in the regular part of the 
operator product expansion, do not give important corrections. In the simplest case where 
we consider the perturbation by a single field where the only primary fields occurring in the 
singular part of its operator product expansion are the identity and the field itself, this gives 
for A(n), 

n^ = -^^ + 0{X^), (3.8) 
dn {2-Kn)y T{hY ^ ^ ' ^ ' 

where C is the three-point coupling. If y > 0, this can be integrated to find the effective 
"exact" coupling Aqo in terms of the TCSA coupling \{n) at level n: 

^-= , .tc^(n) ^H= i^ 2lc^ +0{X')- (3-9) 

This is our one-loop prediction for the coupling constant renormalisation. As we see below, 
this can be improved to take into account the level m of the unperturbed state which leads 
to the replacement of n by n — jn in (|3.9[) . 



4 The ground-state divergence and the energy rescaUng 
4.1 Perturbation theory results 



As in [16 1, the energy rescaling arises as the sub-leading correction to the coupling to the 
identity operator. The bulk case is not as clear-cut as the boundary case, however, as the 
presence of multiple internal channels means that there are small differences in the rescaling 
for states that arise from different representations. We shall also see that the rescaling does 
not necessarily go to zero for all renormalisable perturbations, and for the perturbation by 
a single field, it diverges with n if h > 3/4. We find these corrections by evaluating the 
eigenvalues of the perturbed Hamiltonian to second order. 

We consider the simplest case of the perturbation by a single field ip of weight h with 
coupling A and the scaling operator 

/■27r j/3 

h = {Lo + Lo-^,) + Y,xJ V^{e'')^, (4.1) 
^ Je=o 

where, again, A = (27r)~^A. The eigenvalues of h are the scaling functions and we denote the 
i-th eigenvalue by Cj and take its expansion to be 



e,; 



(A) = ^e,,™A"^, (4.2) 



If /i > 1/2 then one or more of these coefficients will formally be divergent. For example, if 
the unperturbed state \i) is a highest weight state then the first three coefficients are 

Sift = {'^hi - , ej,i = 27rCi^i , (4.3) 



\z\<l \ \Z\ 



= -2vr / — ( i^W{l)v{z)\^) - ) . (4.4) 
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For all states except the the vacuum, \i) = |0), the integral in ej^2 depends in detail on the 
model in question but is divergent if h > 1/2. For the vacuum case with h < 1/2, the third 
and fourth coefficients are given in |^0[ asQ 

= -i(2vr)^2(l - |)7(?/ - 1) , = ^7'(^ " Ihi'i " l)C^^^ , (4.5) 

where 7(x) = r(x)/r(l — x) and A denotes the analytic expression. The expressions ( [4. 51 ) 
can be analytically continued to /i > 1/2 and then agree with the coefficients in the corre- 



sponding TBA calculation. As pointed out in |2C], the TCSA method does not reproduce the 



analytically continued expressions but instead approximates the divergent expression ( [4.4| 
We now demonstrate how we calculate these terms. 



4.2 TCSA results 

The first truncation effects arise in the coefficients €{^2- We shall denote the contribution to 

'[nl 

Ci^rn from the states at TCSA truncation level nhy and the full coefficient in the exact 



TCSA expansion at truncation level n by e"^, so that the TCSA approximation to ej(A) is 

n 

en\) = ^el^~X^, el^ = ^e^l. (4.6) 



k=0 



The term e["J comes from level n intermediate states in the four-point function {i\ip{l)ip{z)\i) . 
In the boundary situation in [16] we could arrange the boundary conditions so that the four 
point function in that calculation was given by a single chiral block; in the bulk this is not 
possible. Instead, we expand the four point function in (4^) over the set of chiral blocks^ as 



(4.7) 





if if) 


{i\^{lMz)\i) =Y,{C^^J? 




3 




j 


' 1 



The contribution from the states at level n in the j intermediate channel comes from the 
coefficient of As in |l§, we find the leading n-dependence of this coefficient by 

expanding the conformal block in powers of (1 — z). 



\ ^ 





] 


[ z 



k 



^(1) 



k 



ip{z) 



^(1). 



— H^iz) 

^ +... 

— i 



(4.8) 
(4.9) 



where Fj^ = -FjA.-i'^f ] are the crossing matrix elements and 



V^(l) 



V'(l) 



ip{z) 



ip{z) 



(1-z) 



~2h 



1+ 2^(1 _ 2)^+0(1 -z) 



(1-Z)-'^(1 + |(1-Z)+0(1-Z)2) 



(4.10) 



■^Note that these differ by powers of {2n)^ as uses the expansion parameter A, not A. 
^See H for details 
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We have assumed that the identity and are the most singular fields in the operator product 
Lp * {p\ if not, there will be correspondingly more terms in ( [4. 91) . Hence we find the leading 
terms arising at truncation level n in the chiral block are 



J I »^ 



1 



^n—hi—h+hj 



Z 



F^^T{n-hi+hj+h) 



+ 



F;r{n~hi+h,] 



T{2h)T{n-hi+hj-h+l) T{2h)r{n-hi+hj-h+r 
This gives the second order truncation level n contribution to the energy as 



+ 



(4.11) 



-i,2 



2{n-hi+hj)V{n-hi-h+hj+lf 



i T{n-hi+h+hj) . T{n-hi+hj) 

-Tjl X 1- i'^.n ^^7T^ H 



r(2/i) 



T{h) 



2r(2/i)2 



2r(2/i)2 



T{h)T{2h) 2T{h) 
Using the crossing properties of the full correlation functions, we find 



4h-5 „3/i-4 „2/^-4^ 



i i i 

so that the second order truncation level n contributions ( [4.12|) are 

-27r 



.(4.12) 



(4.13) 



-1,2 



where 



2r(2/i)^ 



-0,2 



-27r 



2r(2/i)2 ' 2r(/i) 



2V{2hY 



+ (4/i - 3) 



2V{2hf 
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+ 0(n 



E ^/ij , ao = h , 



4/1-5 „3fe-4 „2/i-4 



,71 



n 



lip 



5l 



(4.14) 



(4.15) 



where the second expression in the modular S-matrix Sij holds for a diagonal modular in- 
variant |17]. If \i) is a highest weight state then we can then sum these contributions to 
get 



-0,2 



-i,2 



4,2 - 27r 



n 



Ah-2 



+ 



n 



^^-^ 2/1+1 



4(2/i - l)r(2/i)2 r(2/i)2 4 



+ . . . 



(4.16) 



n 



4h-2 



+ 



n 



4h-3 



2(ctj /ij)+l ^ CiiipCipipipTi 



2h-2 



+ ... 



4(2/i - l)r(2/i)2 r(2/i)2 4 2(2/1 - 2)r(/i)2 

If the TCSA coefficients are converg ent, then (|4.16D gives the leading corrections; if they are 
divergent, they give the leading divergences. We first consider this result for the ground state 
scaling function. 
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For h > 1/2, the leading term(s) in ( [4.16D are divergent. To get a proper estimate of 
the ground state energy from TCSA we need exphcitly to subtract these. We do this for 
the divergent case considered in 0] in section 5.1. As wih be seen, this is very effective at 
producing a good estimate from the TCSA. 

For excited states, we shall only consider the scaled energy gap, = e^ — eo = X] ^i,k^^ and 
its TCSA approximation e" = ^ e^^A^ . Using ( [4.16| ), we can write the TCSA approximation 
as 

n^-'^y h- — 8- n^y 
er(A) = {2hi) + (27r)AC,^, + ( ef^ + 2tt——^^—^ + 27rC,i^a 



"^'^ V ' r(2/i)2 2 ^ '^^'^2yr(/i)2 

= rj,(A) [<5.(Aj„(A)) - 2<Si] + 24^ + 0(\^) , (4.17) 

where 6i = ai — hi. which defines the one- loop energy rescaling rjj and coupling constant 
renormalisation 5n(A), 

The coupling constant renormalisation agrees perfectly with the result of the previous section, 
(|3.9|) . The energy rescaling is a new prediction. It differs from the boundary case as the energy 
rescaling is not exactly the same for each state, because there is at this order a small shift 5i 
which differs between states arising in different representations. Since this is a small overall 
constant shift in each representation, it is only important for the lowest lying states in each 
representation where it does make a noticeable difference, as we see in the plots in section 
p. 

These results also hold, suitably adjusted, for excited states. If \ip) is a state at level m 
in the representation i then the contribution in ( [4.7| ) from the intermediate states at level n 
comes from the coefficient of and following the changes through the net result 

is to replace n in ( [4.18 ) by (n — m), so that for the state l^/)) the rescaling and renormalisation 
are 

''"^^^ " ^ 4r(2/i)2(27r(n-m))2y-i) ' ^"^^^ " ^ 2yr(/i)2(27r(n - m))f ) " ^^'"^^^ 

This alteration is a sub-leading effect in n but is appreciable for the cases we consider when 
n is sometimes quite small. 



5 Tests of the TCSA in the tri-critical Ising model 

One of the first paper to use the TCSA extensively was [Q]. This investigated perturbations 
of the tri-critical Ising model with mixed success. With hindsight, it is easy to see now why 
they obtained good results for the perturbations by the fields of weight 1/10 and 3/80 (with 
fast convergence), mixed results for the field of weight 7/16 (with slow convergence) and very 
poor results for the massless perturbation by the field of weight 3/5. We reconsider this final 
case in the light of our results in section ^.l| . Klassen and Melzer also used the TCSA to test 
predictions of the IR scattering description of the massive perturbation by the field of weight 
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3/5 in [0]. They too found rather poor agreement with TCSA - we are able to improve on it 



greatly using the renormalisation and re-scalings, which we show in section |5.3 . 

Following [^], we denote the massless perturbation of the tri-critical Ising model by 
MA''^^ and the massive perturbation by AiA^^ \ These have been investigated extensively 



and TBA systems found for the ground state(s) in [^2| and [^] respectively. In both cases 
once can define a mass-scale m and a dimensionless system size r = mR satisfying 




For A4A^^ m is in the mass of the kink in this model; for AiA^^' it is an inherent mass-scale 
in the problem and can be identified with the kink mass by analytic continuation in A. 



10.83... |A| 



(+) 



5/4 



(5.1) 



5.1 The ground state energy in A^A 



(+) 



In figure |l(a)| we show the TCSA estimates of the ground state energy (27r/r)eo(r) for the 
massless perturbation AiA^^^ plotted against r for truncation levels 3, 4 and 5. As Lassig et 
al observed, the ground state energy does not appear to be converging to the expected linear 
behaviour of a massless theory. This was explained in as a consequence of the divergence 
of the ground state energy. 



In figure 1(b) we show the same TCSA data but with the leading divergence and sub- 
leading correction given in ( 4.16| ) subtracted, that is, we plot 



^) (e^(A) - i2.r'yx' 



n 



ih-2 



+ 



n 



^^-^ 2/1+1 



4(2/i - l)r(2/i)2 r(2/i)2 4 



(5.2) 



against r, where h 
which is 



3/5 and r is given in (|5.1| 



TT r 
12^ 4 



We also plot the expected IR behaviour 

(5.3) 



where the linear term can be deduced from the perturbative expansion of the TBA solution 
given in 122]. 



As is clearly seen in figure |^, after subtraction of the leading divergence, the TCSA data 
appears to be converging towards the expected IR behaviour. This convergence is dramatically 



improved by incorporation of the leading coupling constant renormalisation as in figure |l(c 
Since this is the ground-state, there is no energy rescaling to include. 



5.2 The energy gaps in M.A 



(+) 

4 



The scaling function gaps in M.A^^^ were also investigated in 0. They found poor conver- 
gence of the first even excitation, but good evidence of the IR fixed point being the Ising 
model. We report here that in fact we have found rather different convergence of the first 
even gap to that presented in 0], something we find hard to explain. 



In figure 2(a) we show the bare TCSA data, which is to be compared with figure 13(b) 
of Ip. As can be seen, the bare TCSA data appears to be converging slowly to its IR limit 
e = 1 for A < 2 but diverging for A > 2; this is in accord with the idea of a fixed point at 



a finite value of A which can be read off from ( 4.18 ) and convergence of the raw TCSA data 
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(a) The bare data 



(b) The subtracted data 



(c) The subtracted, renormahsed 
data 



Figure 1: The ground state energy for the model A4A\^ ' at truncation levels 3 (dotted), 4 
(dashed), 5 (solid) and 9 (dot-dashed) plotted against r together with the leading exact IR 
behaviour (thin solid line). 



only below this fixed point. We also note that our TCSA data is certainly changing more 
slowly with n than was found in 0. As before, the convergence is dramatically improved 



by including the 1-loop rescaling and renormalisation as shown in figure 2(b) . The putative 
fixed point is pushed far off to the right by the renormalisation and the rescaling brings the 
data down below its IR value. 



(a) The bare TCSA data 




(b) The renormahsed, rescaled 
TCSA data 



Figure 2: The first gap in the even sector of AiA^ 
5 (solid) and 9 (dot-dashed), plotted against A. 



(+) 



at truncation levels 3 (dotted), 4 (dashed), 



We would also like to check that the TCSA method gives good results even when the 
renormalisation and re-scalings have become large. To this end, in figure |3| we also include 
a plot of the normalised energy gaps in the (r, *) sectors at truncation level 9. We have 
normalised these relative to the first gap in th (1,*) sector. Assuming that the IR values of 
the gaps are reached around log(A) = 2 and that the first gap in the (1, *) sector tends to the 
scaling value 4 consistent with the identification of the state in the IR as L_2L-2|0), then the 
rescaling function r9(A*) ~ 2.8, which is definitely outside the perturbative regime, They also 
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appear to show that the Ising fixed point is actually reached at about log A ~ 2, in agreement 
with the prediction from the one-loop calculation of a fixed point at a finite positive A. For 
n = 9, the 1-loop prediction for the position of the fixed point is log(A) = 2.45 . 

Despite the size of the rescaling (and the presumably infinite coupling renormalisation if 
the IR fixed point is actually reached) these figures they show good qualitative evidence for 
a flow to the Ising model, the (r, *) sector flowing to the representation. 



(a) The first 19 gaps in the (b) The first 15 gaps in the (c) The first 12 gaps in the 

(1, *) sector (2, *) sector (3, *) sector 

Figure 3: Normalised energy gaps for ^AA^^^ plotted against log(A) at truncation level 9. 



5.3 The energy gaps in the massive perturbation M.A[ K 

The TCSA was also used in in an attempt to verify the description of the IR limit of the 
massive perturbation A4A^^ ^ by a kink model and the approximation of the energy levels 
using the Bethe-Yang equations (we refer the interested reader to Q for details). The Bethe- 
Yang (or BY) equations give approximate finite-size energy gaps for massive models and agree 
with the exact TBA results up to corrections which are exponentially suppressed for large r. 



In figure 4(a) , we reproduce figure 7 from Q showing the scaling functions plotted against 
r. The solid lines the bare TCSA data and the various dotted and dashed lines are the BY 
approximations to these energy levels. There is a general agreement, but we do not think it is 
good enough to confirm the BY energy levels as correct. As Klassen and Melzer say, the main 
problem is not with the BY levels but with truncation errors in the TCSA. In figure 4(b), we 
correct the TCSA data by the 1-loop renormalisation and rescaling formulae we have found. 

There is a dramatic increase in agreement with the BY results for the renormalised, re- 
scaled gaps compared to the bare gaps; they still disagree for small r as is to be expected for 
the approximate BY solutions, but the agreement for middling ranges of r, say 2 < r < 4, is 
excellent. In this range, the exponential corrections to the BY solutions have been suppressed 
and the 2-loop corrections to the renormalisation and rescaling formulae are still small. For 
r > 4 the 2-loop corrections to the renormalisation and re-scalings are large enough to show 
a qualitative difference between the TCSA and BY gaps. 
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(a) The bare data 



(b) The renormaUsed, rescaled 
data 



(c) The fitted data 



Figure 4: The gaps for the massive perturbation AiA\ at truncation levels 5 for figures (a) 
and (b) and 9 for figure (c). In all cases the approximate Bethe-Yang two- and four-particle 
energies are also given as dotted and dashed lines, as calculated in |^]. 



Another noticeable effect of the rescaling is to improve the behaviour of the ground states 
in each sector - the inclusion of the shifts 5i makes a qualitative difference in the convergence 
of the ground states; we show these regions in figure ^ where we plot the energy gaps from 
the bare TCSA data and the RG improved TCSA data together with the difference of the 



ground states in the (1,*) and (2,*) sector calculated using the TBA system of [2C]. There 
is a clear difference between the ground states at r = 5 whereas in figure ^(b) this has been 
substantially decreased and the region of good agreement with the TBA extended from the 
region r < 0.6 to r < 2. 





(a) The bare data 



(b) The renormahsed, rescaled data 



Figure 5: The scaling gaps for the massive perturbation AiA^^ ^ at truncation levels 5 together 
with the exact TBA results for the first gap (dotted). 



We can show that it is only the form of the renormalisation and re-scalings used in figure 



4(b) that are wrong by choosing to fit the TCSA data to the BY data, and so deduce effective 
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re-scalings and renormalisations. We choose to fit the first two gaps in the even sector to 
find these numerical renormalisation and re-scaling, and use these for the rest of the data. 
This empirically renormalised and rescaled TCSA data is shown in figure 4(c) along with the 
BY data, where we have chosen to remove the exponential corrections to the ground states 
in each sector and set them to be zero. The agreement is impressive, even for r as large 
as 12, confirming that the TCSA data is indeed very accurate once the renormalisation and 
rescaling has been taken into account. To give some idea of the size, at r = 12 the empirically 
calculated rescaling and renormalisation functions at level 9 are rg = 0.76 and gg = 0.71. 

We can use the empirically calculated rescaling to check the scaling form predicted from 
( |4l^ ), that is 



r„(A) = r(An(i-2^)/2) , ri.ioop(x) = (l + 



4r(2/i)2(2^)2s/-i 



(5.4) 



In figure ^, Correspondingly, we plot the estimate of and 1-loop approximation to r{x) in figure 
^. We see that there is actually good agreement with both the scaling form and numerical 
value of the function. 




Figure 6: The function r(x) vs. log(x) from the fitting of the odd sector to the BY lines 
at truncation levels 3 (dotted), 6 (dashed) and 9 (dot-dashed) as well as the perturbative 
prediction ri_ioop (solid) 

This can not be said of the renormalisation function gn{x). The inaccuracy of the BY 
approximation for small A or r means that we cannot effectively use it to check the form 
of the perturbative renormalisation function gn{x). For large values of r we see a definite 
disagreement with both the numerical value of and the scaling form of the 1-loop prediction. 
For large values of r, we see gn{x) ~ g{xn~^''^), the same form as the rescaling function. Since 
this is a larger exponent than that predicted from 1-loop, which would be g{xn~^'^), it is 
quite possible that the larger exponent is correct, arising from higher loop corrections, just 
as the exponent for the rescaling comes from a second order effect. The best hope we have 
of checking this in detail is for the exact TBA excited state spectrum to be calculated on the 
cylinder which would allow comparison with ^i.ioop down to small values of r. 



6 Divergence for /i > 3/4 

One of the main predictions of equation ( |4.18| ) for the rescaling functions is the exponent 
of n^^~^ = n^~^y. This goes to zero for large n for h < 3/4 but diverges for h > 3/4, so 
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the rescaling function will diverge with increasing n for h > 3/4. Since this is the leading 
n-behaviour in the scaling function gaps, it suggests strongly that the TCSA scaling function 
gaps will only converge to the exact answer for h < 3/4. This means that any application of 
the bulk TCSA on the cylinder for /i > 3/4 will only be able to predict the ratio of energy 
gaps. As a demonstration, we give here the first few scaling function gaps for MA[^^ with 
h = 3/5, Ah - 3 = -3/5 and for MA^'^'' with h = 4/5, 4/i - 3 = 1/5. As can be seen, 
the gaps decrease with increasing level for A4A^^\ in accordance with the prediction that 
the rescaling function decreases with increasing n. On the other hand, the gaps increase for 
AiA^g~^ , showing no sign of convergence, in accordance with the divergence of the perturbative 
rescaling function. 




Figure 7: The bare TCSA gaps for the (1, *) sectors of (a) MA''^^^ at truncation levels 3 (dotted), 6 
(dashed) and 10 (solid), and (b) A^Ag^^ at levels 5 (dotted), 6 (dashed) and 10 (solid); (e) shows the 
ratios of gaps in Alg^'' at level 10. 



This does not mean that the TCSA contains no useful data - there is good evidence that 
the ratio of the energy gaps still remains physical, even if the gaps themselves diverge. In 
figure \l{c) we show the ratio of the energy gaps 2ei / 63 for the same range as 7(b) , showing 
that they are converging nicely to the expected (*, 1) sector of Mg^g (the dashed lines on the 
right of that plot). 



7 Complications for bulk perturbations on the strip 

So far we have not discussed the TCSA method applied to bulk perturbations on the strip. 
This was first used in |11] where it proved very accurate for the Lee- Yang model (for which 
y = 12/5). We would like to make a similar analysis to that here to find the leading truncation 
effects on the spectrum, but the form of the correlation functions on the strip make this much 
harder and so far we have not been able to find even the leading term such as a coupling 
constant renormalisation. There is however, numerical evidence to suggest that there are in 
fact no such simple renormalisation and rescaling effects on the strip. In figure ^ we show 
the effect of changing truncation level on the scaled energy gaps in the massive perturbation 
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A^A4 ^ on the cylinder and the strip. Since we look at the ratio of energy gaps, the ground 
state energy and rescaling terms are removed, and we should see simply the effect of the 
coupling constant renormalisation. In figure 3(a) this is exactly what we see for the model 
on the cylinder: as the level is increased, the lines all move to the left, suggesting that a 
renormalisation of the coupling is required. In figure p(b) we perform the 1-loop coupling 
constant renormalisation from 4.18| and indeed the TCSA data from the different levels are 
renormalised onto a single set of lines. 

In figure 8(c) , however, we show the scaled energy gaps for the same model on the strip 
with (1, 1) boundary conditions on each edge. In this case, as the truncation level is altered 
there is no consistent movement of the lines to the left or the right, indeed the second gap 
appears to be almost invariant under the change of the level. It is clear that no single coupling 
constant renormalisation will make the TCSA data from the various levels map into a single 
set of lines; the leading truncation effect on the TCSA on the strip is not a simple coupling 
constant renormalisation. 




(a) The normalised gaps for the (b) The normahsed gaps for (c) The normahsed gaps for the 
(1, *) sector of A4^4 ' on the the (1,*) sector of MA^^~'' on even sector of MA^^^ on the 
cyhnder the cyhnder with renormahsed strip with fixed boundary con- 
coupling constant ditions 



Figure 8: The model A4A\ : (a) the normalised gaps plotted vs. r for truncation levels 6 
(dotted), 8 (dashed) and 9 (solid) and (b) the same gaps with renormalised r; (c) shows the 
model on the strip for truncation levels 10, 12 and 16 respectively. 



8 Conclusions 

We believe we have presented good evidence that the leading corrections to the TCSA method 
on a cylinder are a possible ground-state divergence, a coupling constant renormalisation and 
a representation-dependent energy re-scaling. Application of these to the early tests of TCSA 
in and Q have greatly improved the accuracy and reliability of the TCSA method which 
appeared to be very poor in those cases. One by-product has been the important result 
that the TCSA does not actually converge for perturbations with h > 3/4 where there is a 
divergent energy rescaling, but we have presented evidence that the ratio of energy gaps does 
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still converge in this case. 

We would like to point out that the perturbative renormalisation and rescaling use no 
more information than is usually available when using the TCSA - the scaling dimensions, 
three-point couplings and modular S-matrix. In this way we think of them as an improvement 
of the TCSA. They do not require knowledge of the four-point functions or conformal blocks 
which can be difficult to calculate, even for minimal models 

It has also been suggested |24] that the TCSA may have a finite radius of convergence for 



any n. The apparent lack of convergence in n for the first even gap of MA^^'' in figure 2(a) 



might be a sign of this, although since this occurs when the scaling function is greater than 
its IR value, it is also possible as we suggest that this divergence occurs after the IR fixed 
point is reached, which we expect to be at a finite value of n tending to infinity as n increases. 

The results have found now allow one to apply the TCSA in a wide variety of cases where it 
was not thought useful before, in particular we hope to use it to study irrelevant perturbations 
in the light of the results on irrelevant boundary perturbations in |lf 
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